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Abstract: Bosonization is normally thought of as a purely two-dimensional phenomenon, 
and generic field theories with fermions in D > 2 are not expected be describable by 
local bosonic actions, except in some special cases. We point out that 3D SU{N) gauge 
theories on R 1 ' 1 x S\ with adjoint fermions can be bosonized in the large N limit. The 
• key feature of such theories is that they enjoy large N volume independence for arbitrary 

circle size L. A consequence of this is a large N equivalence between these 3D gauge 
theories and certain 2D gauge theories, which matches a set of correlation functions in the 
3D theories to corresponding observables in the 2D theories. As an example, we focus on 
a 3D SU(N) gauge theory with one flavor of adjoint Majorana fermions and derive the 
large-A equivalent 2D gauge theory. The extra dimension is encoded in the color degrees 
of freedom of the 2D theory. We then apply the technique of non-Abelian bosonization 
to the 2D theory to obtain an equivalent local theory written purely in terms of bosonic 
variables. Hence the bosonized version of the large N three-dimensional theory turns out 
to live in two dimensions. 
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1 Introduction 

The path integral representation of the generating functional for a D-dimensional quantum 
field theory with fermionic fields ip, bosonic fields B, and sources J 

Z[J]= j VBV^V^ e ifd D x(jF(B,J)ifi +l C(B,J)) 

can always be manipulated into a form which only has an explicit dependence on the 
bosonic fields B. One simply integrates out the fermions giving 

Z[J] = J VB det(F[B, J]) e lS ^ = J VB e itrio g F(B,j) + iS(B,j) = J VB e iS B (B,j) ^ 

which is a description of the theory in terms of a purely bosonic action Sb(B, J). However, 
the resulting action is in general highly non-local, and this is not what is usually meant by 
the term 'bosonization'. A bosonized representation of a field theory is one in which a set 
of new bosonic variables B' are introduced in such a way that the generating functional 
takes the form 

Z[J}= f VBVB'e iS B( B ' B '^. 
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with a local action S'. Such bosonized representations are known to be available for field 
theories living in two dimensions [1-6] , and are often very useful as they provide a comple- 
mentary description of the physics. As just one example, bosonization makes it simple to 
compute string tensions in a number of 2D gauge theories [1, 7]. 

Unfortunately, above two dimensions local bosonized representations of field theories 
with fermions are not expected to exist in general. In some cases it is possible to approx- 
imately describe the low energy physics of a field theory containing fermion fields using 
local bosonic actions, as is done in e.g. chiral perturbation theory for 4D QCD, or for 3D 
field theories with heavy fermions in [8-15]. While for entirely generic higher-dimensional 
field theories it seems clear that such low-energy bosonization is the most one can hope 
for, it is natural to ask whether there might be special classes of higher-dimensional field 
theories for which exact bosonization is possible in the same way as in 2D theories. 

Our goal here is to point out that 3D SU(N C ) gauge theories with Nf flavors of adjoint 
fermions are an example of such a class of theories, and can be bosonized in the large 
N c limit. In fact, while our discussion will focus on the 3D example for simplicity, in 
the conclusions we will argue that our approach should generalize to 4D gauge theories 
as well. To find the bosonized description of the 3D gauge theory, we move in two steps. 
First, we argue that in the large N c limit, a broad subset of the gauge-invariant correlation 
functions of the 3D gauge theory can be computed in a certain 2D gauge theory thanks to 
the phenomenon of large N c orbifold equivalence. The second step consists of applying the 
technique of non-Abelian bosonization to the 2D theory, and yields a local two-dimensional 
bosonic theory which is large- A^-equivalent to the original 3D gauge theory with adjoint 
fermions. 

The first step in the argument takes advantage of a remarkable property of large 
N c gauge theories known in various contexts as volume independence or Eguchi-Kawai 
reduction [16-18]. For our purposes, large N c volume independence is the assertion that 
for a D-dimensional large N c gauge theory compactified onafl- 2-dimensional spatial 
torus, so that the theory lives on R 1 - 1 x {S^)® 2 with periodic boundary conditions for the 
fermions, there is a non-trivial set of correlation functions that are independent of the torus 
size L, up to corrections that are suppressed by powers of 1/N C , so long as certain symmetry 
realization conditions are met [18, 19]. In particular, so long as volume independence holds, 
correlation functions of a theory on a torus of size L coincide with the correlation functions 
in the decompactification limit L — >■ oo. 

For pure Yang-Mills theories, spontaneous breaking of center symmetry invalidates 
large N c volume independence for torus sizes smaller than a critical size L c set by the 
strong scale of the theory. As a result, the small-volume physics is volume-dependent 
and very different from the physics seen for L > L c , where the large N c theory becomes 
L-independent. However, it was recently realized that the situation is very different for 
theories with adjoint fermions with periodic boundary conditions on the torus, because 
the adjoint fermions prevent the center-symmetry-breaking phase transition from taking 
place. As a result, YM theories with adjoint fermions enjoy the property of large N c volume 
independence for arbitrarily small L [18]. 

The existence of general recipes for bosonization of 2D theories and the phenomenon 
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of large N c volume independence for arbitrarily small L in the theories we consider sug- 
gests a natural strategy for obtaining bosonized descriptions of SU (N c ) YM theories with 
adjoint fermions. We send the spatial D — 2-dimensional torus size to zero, and perform a 
dimensional reduction in a way consistent with large N c volume independence, obtaining 
a two-dimensional field theory which at large N c is equivalent (in a precise sense to be 
defined below) to the original D-dimensional field theory. The application of non-Abelian 
bosonization [3] to the 2D theory then produces a local bosonic theory which has correla- 
tion functions that coincide with those of the original 3D fermionic theory in the large iV c 
limit. 

We hasten to emphasize that the dimensional reduction procedure relevant for such 
theories is not the naive one where one simply truncates the theory to the lowest Kaluza- 
Klein modes, which would of course have no chance of producing a theory equivalent to 
the original large- L theory. Instead, as we explain in Section 2.1, the systematic way to 
perform dimensional reduction for such theories uses the machinery of orbifold projections 
[18], involves the use of gauge- invariant lattice regularization for modes on S 1 , and results 
in 2D theories with unitary scalar fields rather than the Hermitian fields which would arise 
from naive dimensional reductions. 

In this paper we focus on the simplest example of these ideas, which involves studying 
3D SU(N C ) gauge theory with one flavor of adjoint Majorana fermions with bare quark 
mass m on Mr> x S^, which is described by the action 

S 3D = J d 3 xTr (-±^F ap F a P + M] i?j , (1.1) 

where g 3 is the Yang-Mills coupling, ip is a Majorana fermion field and D a = d a — i[A a , ■] 
with A a = A l a ti the SU(N C ) gauge field with field strength F a p, U,i = 1, . . . ,N% — 1 are 
the Hermitian generators of SU(N C ) obeying Trfjij = 2<%> ol, f5 = 0, 1,2, and we use the 
mostly-minus metric convention throughout. In the limit M — > 0, this theory becomes 
J\f = 1 super- Yang-Mills theory. The bosonized theory is obtained from the 3D theory on 
a discretized S 1 direction with a T-site lattice spacing a, L = Ta. We show that in the 
large N c limit, the physics described by Equation (2.1) can also be described by the action 

S= J d 2 x^r(^jl v + ^\D^\ 2 + ^\D^ +m 2 Tr 5 Tr 5 t_^Tr(^)Tr( ff t0t ) | 

+ Nf(g,A), (1.2) 

where [i,v = 0,1, g encodes the bosonized representation of the fermions and is matrix- 
valued field living in the group SU(N c )/Xn c , with TL^ c being the center subgroup of 
SU(N C ), (j) £ SU(N C ) and is related to the component of the original gauge field in the S 1 
direction, in is a mass parameter related to the original quark mass M, and the dimension- 
less parameter c is discussed in Section 2 and Section 5 and has to do with the discretization 
we adopt. Finally, the second line contains an appropriately gauged Wess-Zumino-Witten 
term. 

We hope to make this paper mostly self-contained, and hence organize the exposition 
as follows. In Section 2 we explain the orbifold projections relating the original 3D theory 
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and its large iV c -equivalent 2D form, using the technology developed in [18, 20-23]. On a 
first encounter the phenomena of large N c volume independence and large N c equivalences 
between 2D and 3D gauge theories look quite counterintuitive, and so in Section 3 we work 
through a calculation of the vacuum energy in the 2D theory in perturbation theory, which 
crisply illustrates the way in which the extra dimension is encoded in the color degrees of 
freedom of the 2D theory. The calculation serves as a perturbative verification that the 
vacuum energy of the 2D theory vanishes so long as the fermions are massless, as it must, 
since at m = the 3D theory has J\f = 1 supersymmetry and must have vanishing vacuum 
energy. 

Next, in Section 4 we briefly review the derivation of the technique of non-Abelian 
bosonization of adjoint fermions, focusing on the points which present particular subtleties 
for our application. Finally, in Section 5 we use the previous results to bosonize the 2D 
theory which is large- N c equivalent to 3D YM with adjoint fermions, obtaining a theory 
with the action given in Equation (1.2). In Section 6 we summarize the results and outline 
some promising directions for generalizations and applications. 

Note Added: As this paper was being finalized, a very interesting preprint by Aharony, 
Gur-Ari, and Yacobi appeared proposing an (apparently) unrelated approach to the bosoniza- 
tion of 3D theories of free fermions [24] . 

2 Dimensional reduction of a volume-independent theory 

The starting point of our analysis is establishing the relation between the YM theory with 
adjoint fermions theory living on R 1 ' 1 x and a large- iV c -equivalent theory living on M. > . 
The 3D theory has the action in Equation (2.1): 

S 3D = [ d 3 xTr (-*F a pF af} + $[ip-M]ip). (2.1) 

Jr 1 > 1 xS'1 V *g 3 J 

with a, (3 = 0, 1, 2, with the condition X2 ~ X2 + L and periodic boundary conditions for all 
fields on the S , our convention for gamma matrices in 3D is 7 = 01, 7 1 = 202, 7 2 = io~3, 
and ip is a two-component Majorana spinor. With our conventions, the Majorana condition 
on the spinor components of tp is given by tp = (ipi,iip2) with both tp* = ipi. In the limit 
M — > 0, the action becomes that of N = 1 super- Yang-Mills theory, but it is conjectured 
that the supersymmetry is spontaneously broken unless one turns on a Chern-Simons term 
with a coefficient which is larger than N c /2 [25] . 

Aside from the obvious spacetime symmetries, the theory has an SU (N c ) gauge symme- 
try, and also a global symmetry known as center symmetry. Center symmetry plays a very 
prominent role in the lattice gauge theory and continuum studies of confinement (see e.g. 
[32, 33]), and appears also in the context of string theory, see e.g. [34, 35]. This symmetry 
is somewhat unusual in that its existence as a non-trivial symmetry depends on the topol- 
ogy of the spacetime on which the field theory lives. To recall how such a symmetry arises, 
note that constant gauge transformations that differ by an element of the center subgroup 

For an interesting recent analysis of 3D J\f — 1 SYM theory, see [26], which is based on the approach 
developed in [27-31]. 
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= Ztv,. of SU(N C ) act trivially on all of the fields in the theory, so that one really has an 
SU(N c )/Zn c symmetry acting on the fields [36]. Integrating over the gauge fields amounts 
to taking a quotient by SU(N C ), and the quotient Z = [SU(N c )/Z Nc ]/SU(N c ) = Z Nc is 
the global symmetry, known as the center symmetry. If the theory had lived in M 1 ' 2 , all 
observables in the theory would transform trivially under Z, and one would not really re- 
gard it symmetry worthy of the name. However, on K 1 ' 1 x the situation is different 
because Wilson lines that wrap the S 1 transform non-trivially under the action of Z, which 
acts by changing the periodicity condition on A<i by multiples of elements of the SU (N c ) 
group center. Hence Wilson loops wrapping S 1 transform as 

(TrVe Is i dX2 A ^^) = (n(a/*)> -> w(n(a^)) (2.2) 

where V represents path ordering and oj G ^n c - So Z is a bona-fide global symmetry of 
the theory. 

The theory defined by Equation (2.1) has distinct phases depending on whether the 
ground state is invariant under the center symmetry, with the phases distinguished by the 
expectation values of Trf2 n , |n| < N c . Phases of a theory with unbroken center symmetry 
are known as 'confining', with a linear rising potential between widely separated funda- 
mental color test sources, while phases where the center symmetry is broken are known 
as deconfined phases [32]. Three-dimensional YM with one flavor of adjoint Majorana 
fermions is expected to be confining in large volumes A3L = g^N c L 3> 1, with an unbroken 
center symmetry, which is also known to be true of pure YM theory. The remarkable 
feature of YM theories with adjoint fermions is that they remain confining even at small 
volumes A3L <C 1, as can be seen from a perturbative calculation of the effective potential 
for the eigenvalues of a Wilson loop wrapping the S 1 direction [18]. Such calculations 
show that so long as the fermions have periodic boundary conditions on S , so that one is 
dealing with a spatial compactification as opposed to a thermal one, the effective potential 
is minimized on configurations where the N c eigenvalues of O are evenly spaced on the unit 
circle. This implies 

(Trft n ) =0 for all n / G Z, \n\ < N c (2.3) 

and hence center symmetry is preserved. 

We now explain the construction of the 2D theory which is large- N c equivalent to the 
3D theory. As we will see, it is not given by a conventional dimensional reduction of the 3D 
theory. The most naive dimensional reduction for small L would truncate Equation (2.1) 
to the Kaluza-Klein zero-modes, yielding 



Snaive = L d xTv 



1,1 



■-^f* + \Tv\D^H\ 2 +^(ip + i 7m [H, ■ ] — M)ip 

Z 93 93 



with n, v = 0, 1, 7* = ?7o7i = —72 is the chiral gamma matrix in 2D, and H is the KK 
zero mode of the A2 component of the gauge field. Such a reduction is wrong for 3D YM 
with adjoint fermions, since it amounts to an expansion around a phase with broken center 
symmetry (A2) = 0, and hence (H) = 0. This is not the correct phase for YM with adjoint 
fermions either at large L or small L. 
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A less naive dimensional reduction, which would be consistent with the global sym- 
metry of the 3D theory, would consist of an expansion around a vacuum where A2 (and 
hence H) has a vacuum expectation value (VEV) consistent with center symmetry, for 
instance the Hermitian traceless configuration, (H) = diag (2ir (— N c + 1)/N C , 2tt (— N c + 
3)/N c , . . . , 2ir(N c — 1)/2N C ), leading to a breaking of the SU(N C ) gauge symmetry down to 
U(l) N --K For our purposes, such an approach suffers from two problems, one conceptual 
and one practical. 

The conceptual problem is that expanding around any particular such VEV requires 
gauge fixing, since the values of the entries of the matrix A2 are gauge-dependent. Gauge 
fixing would be an inauspicious way to start an analysis which is meant to work beyond 
perturbation theory, thanks to notoriously tricky issues such as Gribov ambiguities. Since 
our aim is to find a dimensional reduction which would make sense outside of perturbation 
theory, this is a strong reason to avoid such an approach. 

There is also a practical problem with doing a dimensional reduction with a non-trivial 
VEV for A2- Once we fix the VEV in Equation (2.3) for the Polyakov loop, the adjoint 
Higgs mechanism gives a theory where the KK zero modes are separated from the non-zero 
modes by energies of order l/(iV c L), not 1/L as one might have naively expected. The 
factor of 1/N C in the spacing of modes carrying KK momenta arises from the spacing of 
the eigenvalues of H demanded by center symmetry, and indeed is the heuristic reason for 
the volume independence of gauge theories in center-symmetric phases. In the 't Hooft 
large N c limit where \%dL is fixed as N c — > 00, which is the most interesting one from 
the perspective of the 3D theory, one would have to keep the whole infinite tower of such 
states, and the resulting '2D' theory would not truly look two-dimensional. 

To construct the desired 2D theory, we need an approach in which a manifestly gauge- 
invariant dimensional reduction can be defined, while allowing for the preservation of center 
symmetry. It turns out that the necessary approach is to first discretize the S 1 direction 
on a lattice with T sites and lattice spacing a, so that L = Ta, and borrow the standard 
technology of lattice gauge theory to keep the resulting theory manifestly gauge invariant. 
Then the volume-independence of the theory translates into the independence of the theory 
on the number of lattice sites T. With such a formulation, there is a very natural and 
simple way to define dimensional reduction to 2D that does not require gauge fixing or a 
commitment to expanding around the wrong vacuum. One simply defines the desired 2D 
theory to be given by discarding all of the lattice sites except for one, which removes all 
link fields except for one, yielding a theory living on R ' [37]. Of course, one must then 
explain why such a prescription would yield a theory which is equivalent at large N c to the 
original 3D theory. The reason for the equivalence is because such a 2D theory is connected 
to the 3D theory by orbifold projections which imply a large N c equivalence [18]. 

Before explaining the orbifold equivalence, we write down the discretized theory. Bor- 
rowing language from the story of dimensional deconstruction [38], the discretized theory 
is defined by encoding A2 >n = ^2(^2 = na) in unitary link fields U n = e mA2 '",n = 1, . . . , T, 
and the Lagrangian takes on form of a quiver gauge theory with T SU(N C ) nodes with 
4> n transforming as (F,F) under [SU(N c ) n -i, SU(N c ) n ]. The fermion discretization is 
more subtle and described below. Our discretization of the 3D theory is defined by the 
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Lagrangian 

Siat = J ^ x a yZ Tr 



— (F { ) F l ^ v A — ID^-I 2 



-^(iV-D* - rn)ipi + ^^i73 jWi+i^f - <^_i^-i^-i} (2.4) 

where = 8^ + - Djrf = + t[Aj,, # 

The fermion terms and the parameters m, r and c deserve explanation. On the second 
line, there is a standard fermion kinetic term for 2D fermions, and a term which becomes 
the kinetic term along X2 in the continuum limit X^d - 0. As is well known, however, 
naive discretizations of fermions suffer from fermion doubling, meaning that with only the 
terms in the second line we would end up with a continuum theory with two flavors of 
fermions, not one. 

We deal with this in the simplest and least sophisticated way, by including a Wilson 
term with coefficient r in the third line. Eventually, we will set r = 1 as is standard in 
the lattice literature. The Wilson term, which in the continuum limit approaches a'tpD^ip, 
serves to break the lattice doubler symmetry [39] which ensures fermion doubling for the 
naive fermion action with r = 0. The Wilson term gives the doubler fermions masses of 
order 1/a, so that they decouple in the continuum limit. However, the Wilson term also 
breaks parity symmetry (the analog of chiral symmetry in 4D), and as a result the physical 
fermions also pick up an additive mass shift of order 1/a thanks to interactions. To get a 
single flavor of fermions with an adjustable mass, we add a fermion mass term m, which 
has to be tuned to produce the desired 'physical' mass, and m / M. 

Another subtlety we must deal with has to do with renormalization of the discretized 
fermion kinetic term, and this is the reason the Wilson kinetic term appears in Equa- 
tion (2.4) multiplied by the dimensionless parameter c. The point is that quantum correc- 
tions will renormalize m, c, and the coefficient of |-D</>| 2 - Above we have already discussed 
the renormalization of m, and the renormalization of the coefficient of |-D</>| 2 can be viewed 
as simply a renormalization of the value of a. To make sure the fermions see the same 'speed 
of light' in the S\ direction as the gauge bosons, however, one must tune c appropriately, as 
is well-known in the literature on anisotropic lattice actions, see e.g. [40-46]. At tree level 
in perturbation theory, the correct value of c to get a Lorentz-invariant continuum limit 
would be c = 1, but in general c 7^ 1 to get the continuum limit we want. Tuning c gives 
us control over the only marginal operator 2 consistent with the symmetries which could 
break Lorentz invariance, and hence will be enough to give a Lorentz-invariant theory in 
the continuum limit 3 . 



2 The fact that one does not need to introduce separate c parameters for each term in the sum over 
lattice sites is due to the discrete lattice translation symmetry. 

3 We are grateful to Francis Bursa, Tom Hammant, and Matt Wingate for very useful discussions about 
these issues. 
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Finally, we note that while in 4D gauge theories the tuning of parameters analogous 
to m and c necessary to get the desired continuum limit is in general non-perturbative, 
the situation is better in 3D gauge theories, since 3D YM theory is super-renormalizable 
[47-53]. Consequently it is generally possible to determine the necessary tuning of m,c 
by finite-order analytically calculations in lattice perturbation theory, and doing so for the 
particular theory we are considering here is an interesting direction for future work. 

2.1 Orbifold projection 

Having explained the discretization of the S 1 direction of the 3D theory, we are finally in 
a position to define the dimensional reduction down to two dimensions and to explain the 
reason for the large iV c equivalence between the 3D and 2D theories. 

The point is that discarding all lattice sites except for one can be viewed as an orbifold 
projection. For our purposes, an orbifold projection is a quotient of a 'mother' theory by 
some discrete symmetry G, which defines a daughter theory. In the case of interest to 
us, it is also possible to define an orbifold projection of the daughter theory that yields 
back the mother theory, but with a reduced number of colors. Remarkably, for theories 
linked by orbifold projections, it has been shown that the planar diagram expansions of the 
theories coincide to all orders in perturbation theory 4 , leading to a large N c equivalence 
between correlation functions which survive the projection [20, 21]. In fact, in [22, 23] it 
was proved that for theories with matter in two-index representations, large N c orbifold 
equivalence holds non-perturbatively provided the symmetries used in these projections are 
not spontaneously broken. 

To apply the orbifold projection/equivalence machinery to our case, note that the 
discretized 3D theory has a discrete translation group T s i = Zp. The orbifold projection 
to the single-site 2D theory is simply a projection based on the Zr translation symmetry 5 . 
Such an orbifold projection discards all of the nodes of the quiver except for one, leaving 
only one link field, which becomes an adjoint unitary scalar. The action of the 2D daughter 
theory obtained from the orbifold projection takes the form 

S 2D = [ (fxTrl-^F^F^ + -^\D^\ 2 + ^{i]p-m)^+ 

-^Ah*{i>A ] } + ^A][^,^]) (2.5) 

where g 2 is the two-dimensional gauge coupling related to g% via g\ = g^/a, 4> is a uni- 
tary scalar living in SU(N C ) and is the remaining link field, a = 0, 1, 7* = 27071 = 0-3, 
the fermion fields have been rescaled by a factor of a to make sure they are canonically 
normalized in 2D, and we have included the rescaling of the projected action by a factor 
of 1/r necessary for the 2D theory to have the same correlation functions as the original 
3D theory in the large N c limit. The first fermionic term is simply the usual kinetic term 

4 Provided that the action of the daughter theory is rescaled by a factor of rank(G) -1 before one tries 
to compare the surviving correlation functions. 

5 Reduction to theories with more than one site would involve orbifold projections based on non-trivial 
subgroups of the discrete translation group 
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for an adjoint Majorana fermion with bare mass m in two dimensions. The second term 
is the remnant of the fermion kinetic term in the S 1 direction of the 3D theory, while the 
third one is the remnant of a Wilson term in the 3D theory. 

The 2D theory defined by eq. (2.5) has a SU(N C ) gauge symmetry, with tjj and <f> 
transforming in the adjoint representation. There is a G c = Zjv c global symmetry that acts 
as 4> — * ^(f) with w G Zp, with the other fields transforming trivially. This is the remnant of 
the center symmetry of the 3D theory 6 . There is also a discrete symmetry acting as either 
ipL —^L or ipR — > —tpR, provided that at the same time one sends <fi — > (fr,m — > —m 
and r — > —r, which will be very useful in Section 5. 

Large N c orbifold equivalence implies that the 2D theory will be equivalent to the orig- 
inal 3D theory only so long as discrete translation symmetry in the 3D theory is unbroken 
and the center symmetry of the 2D theory, which is probed by the order parameters (Tr(/> n ) 
is unbroken [18, 19]. The first condition is not expected to be violated in gauge theories 
without chemical potentials. Given the effective potential calculations in theories with ad- 
joint fermions cited above [18], one might suppose that the second condition is sure to be 
met thanks to the presence of adjoint fermions in the theory. However, the realization of 
center symmetry in the 2D theory is somewhat subtle [37, 54-56]. 

The 2D theory is a gauged non-linear sigma model. Non-linear sigma models without 
gauge fields are renormalizable in two dimensions (and indeed can be exactly solved [57- 
59]), and of course gauging does not break renormalizability. However, the gauging breaks 
the SU(N c )l x SU(N c )r symmetry of the sigma model, which acts as <fi — > Lcj)W , down 
to the diagonal SU(N C ) subgroup L = R. As a result one can add a very large number of 
marginal terms to the Lagrangian in Eq. (2.5) consistent with its remaining symmetries 7 , 
namely 



Clearly, these terms amount to a potential for the eigenvalues of eft, and hence by tuning 
these coefficients one can affect the realization of the center symmetry. Whether or not 
such terms are added to the original action, they will be generated by quantum fluctu- 
ations, which will lead to a renormalization of the values of d n>mi .... If one chooses 
not to add the above "double-trace" potential to the action and lets one be generated just 
from quantum fluctuations, the result will obviously be regularization scheme dependent. 
Amusingly, as shown by [54], a complete discretization of the 2D theory on a lattice with 
the same spacing in all directions and a Wilson fermion action with no tree-level double- 
trace potential produces an effective potential of the form V with coefficients that drive 

6 Provided N c — NT for some integer N, which can certainly be arranged for any desired T if N c is large, 
one can define an orbifold projection based on a Zr subgroup of G c and its embedding into SU(N C ) in such 
a way that the resulting daughter theory is again a discretized 3D theory on V sites, but now with N colors 



V = L 4 c n Tr|fi™| 2 + . . . which is consistent with all of the symmetries of the theory but may affect the 
realization of the center symmetry of the theory. These are simply the double-trace deformations introduced 




(2.6) 




in [33, 60, 61]. 
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the theory to a center-symmetric vacuum. If one uses some other regularization that does 
not share this fortunate feature, one could still tune the coefficients in V to arrange for the 
center symmetry to be preserved [56]. Hence large N c orbifold equivalence can be arranged 
to be valid for the 2D theory we will be studying, and in what follows we will assume that 
this has been done. 

The claim that Equation (2.1) and Equation (2.5), possibly deformed by the addition 
of V, define equivalent theories at large N c may raise a number of natural questions, which 
we should address before discussing bosonization. 

1. If the 2D theory is supposed to be equivalent to the 3D theory at large N c , where is 
the extra dimension encoded in the 2D theory? 

2. Does the 2D theory know about the fact that in the continuum limit the 3D theory 
has M = 1 supersymmetry if m = 0, with powerful consequences for e.g. the vacuum 
energy of the theory? 

3. In the original 3D theory, at least for large L, there is certainly no need to add 
any potential of the form of V to the action. But if the 2D and 3D theories are 
to be equivalent, and a non-trivial V needs to be included in the definition of the 
2D theory, it better be the case that the common-sector physics in the 2D theory 
does not depend on the values of the parameters in V so long the theory is in a 
center-preserving phase. But how does that work? 

In the following section we will investigate the second point by a study of the spectrum 
of our theory at the leading order in perturbation theory around the center-symmetric 
vacuum. This will allow to see the fate of supersymmetry by calculating the vacuum 
energy in the 2D theory in perturbation theory. As it happens, the calculation also serves 
to illustrate the answers to the other two questions in a particularly simple context. (Of 
course, the answers to the first and third questions have already been addressed in the 
seminal works on volume independence [16, 18, 33].) 

3 Spectrum of modes and the vacuum energy 
3.1 Perturbation theory and the spectrum 

To illustrate how the dimensional reduction works, in this section we explore its properties 
in the simplest possible context. Specifically, we will compute the spectrum of our 2D 
theory in perturbation theory, and compare to expectations from 3D. The technical goal 
of the analysis is to verify that in an appropriate limit of the 2D theory corresponding to 
the emergence of M = 1 supersymmetry in the 3D theory, the vacuum energy vanishes. 

An immediate problem we must face before beginning the calculation is that as far as 
observables like the spectrum is concerned, in the large N c limit with the 't Hooft coupling 
A3 = <?3 and circle size L held fixed (more precisely holding the dimensionless quantity 
A3L fixed), the theory is necessarily in a strongly interacting regime. After all, as long as 
continuum volume independence holds, which it does in such a limit, the physics at any L 
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will be large N c equivalent to the physics for L — > oo, where we expect the spectrum to be 
inaccessible to weak coupling calculations. 

However, it is possible to find a regime where the large N c 3D theory becomes weakly 
coupled. The trick is to consider a non-'t Hooft large N c limit with L ~ N~ , so that 
X^LN C is fixed as N c — > oo. The downside of such a limit is that it brings us outside the 
regime of validity of continuum volume independence, since the properties of the resulting 
theory will depend on L [62]. The positive aspects of considering such a regime are that 

1. If X^LN C <C 1 one expects perturbation theory and semiclassical approximations to 
become reliable guides to the properties of the theory. 

2. Center symmetry is expected to be preserved for any X^LN C (as can be checked 
explicitly if X^LN C <C 1), and the theory will still be in a confining phase 8 . 

This suggests that if we do a perturbative calculation in the 2D theory and take large 
N c equivalence to the 3D theory for granted, we should only expect the calculation to be 
controlled for X^LN C <C 1 9 . However, once the 3D theory is out of the volume-independent 
regime, it less heuristically obvious that the large N c limit of the 2D theory should continue 
to capture the 3D physics. In particular, the non-perturbative arguments for large iV c 
equivalence presented in [22, 23] were constructed for 't Hooft large N c limits, not the 
large N c limit we are driven to consider. However, the perturbative arguments of [20, 21], 
which match the planar Feynman diagrams between the two theories, will still apply. So 
at least in perturbation theory, the dimensionally reduced 2D theory we have constructed 
will still be large N c equivalent to the 3D lattice theory (and hence the continuum limit 
of the 3D theory as well) even when X^LN C <C 1. Our calculations in this section will 
be perturbative, and so this is good enough for our present purposes. It is tempting to 
speculate that large--/V c equivalence for XsLN c <C 1 will hold non-perturbatively as well, 
but an exploration of this is beyond the scope of the current investigation. 

With this subtlety addressed, we proceed to the calculation of the spectrum in per- 
turbation theory 10 . As in any perturbative calculation, we must start by choosing an 
appropriate gauge and vacuum to expand around. As we have been stressing, the presence 
of adjoint fermions and an appropriate choice of V will guarantee that the reduced theory 
vacuum will preserve the center symmetry. For our perturbative calculation, we choose a 
gauge in which the unitary scalar <f> takes the VEV 



fl ■■■ \ 

w : 

: '•• 

\0 ••• co N ~ l J 



= n. (3.1) 



8 See [62-64] for discussions of how this happens in the parallel 4D case. 

9 For other recent perturbative studies of large N c gauge theories in center-symmetric phases, see [54- 
56, 60, 61, 65-67]. 

10 A pioneering example of such calculations in center-symmetric vacua, see [68]. 
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We will determine the spectrum of the theory to leading order in perturbation theory by 
expanding <p around this vev: 



--n e iH ~ n (1 + iH + 0(H 2 )) 



(3.2) 



where we introduced the Hermitian matrix H, which is assumed to have zero expectation 
value. 

The covariant derivative for <p introduces some mixing terms between and H. 
These mixing terms can be removed by working in i?g gauge, which amounts to adding 
an appropriate gauge-fixing term and ghosts to the Lagrangian. The gauge-fixing function 
and Lagrangian are 



92 



a 2 92 



£ gf = -Tr (G 2 ) . 



(3.3) 



Under infinitesimal gauge transformations parametrized by a = a z ti, the fields appearing 
in G transform as 



5H = - (stfan-a) (1 + iH) - i[a, H] , 5A fl = D fl a, (3.4) 
and from here it is straightforward to compute 5G/5a, obtaining the ghost action 



C 



-Tr 



—cd a D^c--^— (cti)cQ. + cQ.cQ ] - 2cc 



92 



a 2 92 



(3.5) 



To read off the spectrum of the theory to the leading order in perturbation theory 
we simply need the quadratic part of the action. The interaction terms will make small 
contributions so long as X^LN C <C 1. After canonically normalizing all of the fields, the 
quadratic part of the action can be written 



S, 



quad 



d 2 x l A 



x fM nm 
~\~ C-nm 



1 _ I ] d»d v + sin' 



-□ ~ sin 

a 1 

— U k sin 

cr 



+ ic — sin 
a 



7r(n — to) 
N 

7r(n — to) 
N 

2-7r(n — to) 
N 



2 / 7r(n — to) 



9 



A„ 



(3.6) 



H m n~\~ 



2rc . n 

to -| sin 

a 



7r(n — to) 
N 



and since to the order to which we are working the interactions between the modes are 
being ignored, we should set c = 1, which we do for the rest of this section. 

At this stage it becomes clear how the spectrum of the theory can be independent of 
the coefficients in V so long as the vacuum is center symmetric. Consider for instance the 
term 

^Tr(0)Tr(0t) 
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in V, with c\ = 0(N®). Expanding <f> around fi, making use of the fact that H is Hermitian 
and traceless, and switching to canonical normalization, this term makes the contribution 



ci 



N r a 2 



Tr{nH)Tr(&H) - -Tr(fiiT 2 )Tr(fi t ) - -Tr(& H 2 )Tr(n) 



to the quadratic action. If it were the case that f2 = 1at c x7V c , then one would have Trfl = N c , 
and the second and third terms in the expression above would shift the mass of the H 
modes at leading order. However, if Trfi = 0, as is the case so long as center symmetry 
is preserved, the second two terms above vanish, and the first term makes a contribution 
to the mass only at order 0(N^T 1 ) relative to Equation (3.6). It is simple to see that the 
same results hold for the other terms in V: as long as the theory is in a center-symmetric 
vacuum, the values of the coefficients in V do not affect the center-symmetric observables 
in the theory to leading order in the 1/N C expansion. This is a simple illustration of the 
general argument to this regard given in [33]. 

Coming back to our action, Equation (3.6), we can read off the dispersion relations of 
the various modes in the theory. The longitudinal W-bosons (A^ ~ p^) have mass 

„2_„2 , t 4 ^; n 2 ( <n- m 
'a 2 I N 



Po=Pi + t-2^[ V Ar ' , (3-7) 



and the same relation holds both for the scalars H and the ghosts c, c, so they all decouple 
in the unitary gauge £ — > oo. The transverse W-bosons (Ap ~ e^uP u ) are the physical ones: 



/ \2 2 2 2/ /t \' t — / \ i \ 

€ h {pi,n,m) = = + -J sin , (3.8) 



4 . 2 / vr(n - m) 

' -77 sin 

a? \ N 

and the Majorana fermion dispersion relation gets contributions from the Wilson term and 
the bare mass: 



/ n2 _ 2 2 , c 2 ■ 2 ( 27T (n - m) , 
e f (pi,nm) =p =p 1 + ^sm + 



a 2 \ N 

3.2 Vacuum energy 



2rc . 2 / 2it(n — m) x 



m H sm 



a \ N 



. (3.9) 



With the dispersion relations in hand, we are finally in a position to write down an expres- 
sion for the vacuum energy density to leading order in an expansion in \^LN C : 

E vac = E h + E{ (3.10) 

= E /W P ? + W=^)+ (3.1D 

\<n<m J V V / 



Kn<m 



2rc 9 / 2ir(n — m) 
m H sin ( — 



n 2 



where the boson, E b = ^ P:njm eb(p,n,m), and fermion modes, E f = Yl P , n ,m e b(P> n ; m )> 
contribute as usual with opposite signs. We note that the Higgs and longitudinal bosons 
cancel precisely the ghosts contribution to the vacuum energy since the ghosts come as 



-13- 



complex scalars with fermionic statistics. The p\ integrals in the expressions above diverge 
for large p±, and the expression is to be understood as having an implicit UV regulator 
about which we will have more to say shortly, while the sums over n, m are finite. 

Clearly, E vac ^ for generic values of the parameters. How can we reconcile this 
with our 3D expectation that E vac = for massless fermions, on account of the M = 
1 supersymmetry that emerges at m = 0? The key is to recall that the 2D theory is 
supposed to be large- A c equivalent to a 3D theory on a lattice in the circle direction, not 
directly to the continuum R ' x S 1 theory. The lattice discretization we have chosen breaks 
supersymmetry, most directly thanks to the breaking of the U(l) translation symmetry to 
its Zr subgroup. In an appropriate continuum limit, however, we expect the restoration 
of the full translation symmetry, up to corrections proportional to positive powers of the 
lattice spacing. 

Hence to see the fate of E vac , we should see how Equation (3.10) behaves in the contin- 
uum limit. Indeed, we recall that for field theories on a spatial lattice, the dependence of 
the dispersion relations on the momenta is trigonometric, and the momentum integrations 
are restricted to the Brillouin zone. In particular, we recall that for a free scalar field with 
mass M on M 1 ' 2 , with the X2 direction on a lattice with spacing a, the dispersion relation 
is [39] 

^ = M 2 +p 2 + lsin 2 (^) (3.12) 

with p2 € (— 7r/a, n/a]. For small ap\ this is just p$ = M 2 +p\ +p\ , but there are corrections 
in powers of a, and for large momenta (~ 1/a) the the propagator is very different, as one 
could expect since the lattice acts as a UV regulator for p2- 

Armed with this observation, we can see the mapping of parameters necessary to 
interpret Equation (3.10) in terms of the 3D theory. Consider first the contribution to the 
vacuum energy from the gauge-Higgs-ghost sector, which takes the form of 

*= £ l*J* + £**(«zjM. (3.13) 

l<n<m V 

We now suppose that N c = YN for some r,A c ^> 1, as can certainly be arranged for 
sufficiently large N c , and define 

o = ^P2= ]\t " = T ^> 3 - 14 

2 iV c JS c a L N 

where we recalled the relation L = Ta between the lattice size and the lattice spacing. An 
expansion of in ap2 <C 1 now yields 

N c 

E h = Yl dp 1 ^p 2 +p 2 ,[n,m} + .... (3.15) 

l<n<m 

and we finally see that it is the color sums in E^ which encode the third dimension. For 
instance, if we consider n, m such that |n — m\ mod = 0, we pick up the contributions 
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to E\y which would arise from a KK tower for a scalar on a circle of circumference L. 
But of course the fact that the sum also includes momenta with |n — m\ <C N means 
that E\y also includes contributions from momenta much less than the naive KK scale, 
so that the gap between the zero modes and the non-zero modes is actually j^. This 
is precisely the same as what happens in a perturbative calculation of the gauge-Higgs- 
ghost contribution to the vacuum energy the 3D theory in a center symmetric vacuum. In 
both the 2D and 3D theories the KK scale in a center-symmetric vacuum becomes jj^, 
with the momenta becoming continuous and gapless in a large N c , N limit with L fixed, 
which is also the limit in which a perturbation theory becomes unreliable, since it is the 
volume-independent regime. On the other hand, if one takes a non-'t Hooft large N c limit 
with L ~ -/V" 1 , then the KK scale becomes stabilized and order 0(N®), and both theories 
become weakly coupled. The fact that the matching is between theories with different 
numbers of colors N c , N is expected from orbifold equivalence (and falls out naturally from 
volume-expanding orbifold projections), and harmless so long as N C ,N ^> 1, since the 
spectrum must become N, ^-independent so long as N, N c are large enough and the large 
TV" limit exists. 

Finally, we consider Ef, and use the same parameter identifications as in the boson 
case to see the connection to the 3D theory. As is usually the case, the fermion case is 
more subtle. First, to the order to which we are working, we must take c = 1, and will do 
so for the rest of this section. If we also take m = r = 0, then the extra factor of two in 
the the argument of the sin 2 in Ef together with the limits on the color sum implies that 
the near-continuum expansion of the sin 2 includes contributions both from near \ap2\ ~ 
and \ap2\ ~ Hence 

N c 

E { (m = 0,r = 0) = -2 ^ / dpi ^ p\ + p 2 2 [n, m) + . . . . (3.16) 

l<n<m 

and Eb + Ef ^ 0. Indeed, the vacuum energy turns out to be Eb + Ef = —Ef, and it is 
as if there were two flavors of Majorana fermions in the theory, rather than one. In fact, 
this is exactly right. If r = 0, the 3D lattice theory to which the 2D theory is large- iV c 
equivalent suffers from fermion doubling, and for m = there are two physical 'flavors' of 
fermions rather than the one flavor one thought one was discretizing. 

If we turn on r > 0, things change dramatically. The expansion of the trigonometric 
factors in Ef near ap2 ~ vr picks up an order 1/a 2 mass-squared term proportional to r, 
while at the same time the Wilson term makes no contribution to the mass of the modes 
with ap2 ~ to this order in perturbation theory. Hence the doubler mode becomes 
much heavier than the 'continuum' mode we wanted, and so long as the momentum cutoffs 
implicit in E v!)iC are lower than 1/a, the contribution of the doubler mode to the continuum 
physics can be neglected. Of course, this was precisely the motivation for including the 
Wilson term in the lattice theory in the first place: Wilson terms prevent fermion doubling. 

Hence we see that to the order to which we are working and with a sensible choice of 
UV cutoff, we find 



£?vac(m = 0,r ^0) = 0, 



(3.17) 



as one would expect for a theory which has emergent supersymmetry in the continuum 
limit for a massless Majorana fermion. If we were to go to the next order in perturbation 
theory, we expect that for r ^ 0, to get the physical fermion mass to be zero the bare 
fermion mass m would have to be tuned appropriately, since the Wilson term breaks the 
same symmetries as the mass term and hence induces additive mass renormalization [39]. 
At the next order in perturbation theory we would also have to adjust the value of c away 
from c = 1. These observations serve to re-emphasize the points we have made before, 
which are that in general m ^ M, 1. 

4 Bosonization of adjoint fermions 

To bosonize the theory described by Equation (2.5), we need to know the bosonization 
'dictionary' for massive adjoint fermions coupled SU(N C ) gauge fields and scalar fields. A 
bosonization recipe for fermions with non-Abelian color or flavor symmetries was invented 
in [3], and rederived using a constructive path integral framework in [4-6, 69]. These 
path integral derivations of the bosonization dictionary are based on the beautiful explicit 
calculation of the functional determinant of fermions coupled to gauge fields of Polyakov 
and Wiegmann [57, 70]. In Section 4.1, we follow [57, 70] to compute the functional 
determinant for adjoint fermions coupled to SU(N C ) gauge fields. Once this is obtained, 
the methods of [4-6, 69] yield the bosonized representation of adjoint fermions coupled to 
gauge fields without any mass or Yukawa terms. In Section 4.2 we explain how to take 
into account mass and Yukawa terms, which is necessary to enable the bosonization of the 
theory described in Section 2. 

4.1 Evaluation of massless fermion action 

We first suppress the mass terms and scalar couplings, and consider free Majorana fermions 
ip in the adjoint representation of SU(N C ), described by the Lagrangian 

£ adj =Tr[#y Q £> a # (4.1) 

The SU(N C ) gauge fields Aj^ = A'tta, a = 1, . . . , N 2 — 1 in this expression are to be thought 
of as background source fields with local SU (N c ) transformation, and the path integral over 
the fermions defines the generating functional 

Z[Af] = J V^e^ d2xC ^ (4.2) 
= det[C(ij a d a + A^)} 1 / 2 = expi [Tr log W(A S ^)] (4.3) 

where the square root appears because we are dealing with Majorana fermions 11 . Our aim 
is to find an explicit expression for Z[A^ U ] as the exponential of some bosonic local action. 

Before discussing the derivation of the bosonic representation of Z[^j u ] we must dis- 
cuss an important subtlety. The issue is that Z[A^ U ] does not describe all of the correlation 

11 Strictly speaking the result of integrating over Majorana fermions is the Pfaffian, which can differs from 
the square root of the determinant by a non-trivial sign, but this subtlety will not be important for us. 
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functions of currents that could be denned for Equation (4.1). The reason is that in the 
absence of SU(N C ) gauge fields, Equation 4.1 actually has an SO(N 2 — 1) global sym- 
metry under which fermions transform in the vector representation, and the theory has 
a large number of conserved currents that do not couple to the defined above, which 
only parametrizes an SU(N C ) subgroup of SO(N 2 — 1). Hence one could try to bosonize 
adjoint fermions by finding a bosonic representation of ^[^4^°], where ^4^° are background 
SO(N 2 — 1) gauge fields. To obtain the Z[^4^ u ] generating functional above, one would 
then have to set to zero all of the components of A^ that do not live in some chosen 
SU(N C ) subgroup of SO{N? - 1). 

However, given our goal of describing a theory with a bona-fide SU(N C ) gauge symme- 
try, for which we only actually need Z[^4^ u ], the 50(iV ( ? — l)-symmetric approach would be 
overkill. Instead, we will pursue the much more efficient approach of directly determining a 
bosonized description of Z[A^ U ], so that we only explicitly work with a bosonized descrip- 
tion that describes the degrees of freedom that play a role in the SU(N C ) gauge theory of 
interest. The result of this approach will be essentially the same as the one reached in [71]. 

These issues can be understood more explicitly in a much easier setup. If we consider 
a theory of just three free Majorana fermions then at the free level we would have an 
SO (3) l x SO(3)r symmetry, and the bosonization of the theory would give rise to an £0(3) 
WZW model [3]. This amounts to working out an expression for the fermion determinant 
coupled to an 50(3) gauge field, and then using the techniques of [4-6, 69] to get the 
bosonized theory. However, suppose we are actually interested in gauging an 50 (2) y = 
U(l) subgroup of the fermionic theory. On the fermionic side, this is simply a theory 
of a massless Dirac fermion charged under a U(l) gauge field and a neutral Majorana 
fermion; a gauge-invariant mass term could of course be added. In particular since the 
neutral Majorana fermion does not interact with the Dirac fermion or with the gauge field, 
correlation functions involving both the Dirac and Majorana fermion factorize in a trivial 
way. If one wanted to bosonize such a fermionic theory, one could simply use Abelian 
bosonization to capture the correlation functions of the Dirac fermion, and the result is 
a theory of a scalar field coupled to a Sine-Gordon potential and a gauge field. This 
corresponds to working out an explicit expression for the fermion determinant coupled to 
U(l) gauge field, and then applying the techniques of [4-6, 69] to get the bosonized theory. 
Of course, the same physics is contained in the 50(3) WZW model coupled to a U(l) 
gauge field, but the encoding is rather complicated since the 50(3) boson representation 
describes the full set of S0(3)l x SO(3)r correlation functions. But so long as one is 
only interested in the non-trivial gauge-invariant correlation functions involving the Dirac 
fermion field, the Abelian bosonization is equivalent to the full non-Abelian bosonization, 
with the virtue of being far simpler. 

In what follows, then, we seek to obtain an expression for Z(Aj^) directly, following 
the methods of Polyakov and Wiegmann. The method of [57] is based on the construction 
of solutions to the Ward identities of the theory in terms of bosonic variables. To begin, 
note that the classical action has the vector and axial flavor symmetry SU {N)y x SU (N)a, 
with tp transforming in the adjoint representation of both symmetry groups. The vector 
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current is by definition sourced by An'. 



Crucially for what follows, the two-dimensional identity j^j 11 = e^ Vr ^ v implies that the 



axial current is = en U J v,u , which allows us to rewrite relations for J A in terms of J v . 



As is well-known, is a conserved current, while the axial current is anomalous, and the 



currents obey the equations 

d»jX + i[A",jV]=0 (4.5) 
d»J* + i [A», = ^ {d v jv + i [A v , Jj]) = ^F, u (4.6) 

We now follow the method of Polyakov and Wiegmann [57] to find an expression for W{A fjL ). 

The approach is technically simplest in lightcone coordinates x± = xq ± x\, to which 
we now shift. We can rewrite Equation (4.1) in a SU(N c )l x SU (N c ) r manifestly invariant 
way: 

£ adj = Tr[V2 (id+ - [A + , .]) + P R {id- - [A., .]) . (4.7) 

Clearly the vectorial and axial SU(N C ) gauge fields can be recovered by suitable combina- 
tion of left and right SU(N C ) gauge fields. We now rewrite A±, which live in the lie algebra 
of SU(N C ), in terms of group- valued variables: 

A + = -i g^d+g, A_ = -i h^d-h (4.8) 

and g, h live in the group SU(N C ), and pure gauge configurations correspond to the choice 
g = h. In fact, any elements of SU(N C ) differing by multiplication by an element of the 
center G c subgroup of SU(N C ) will lead to the same A±, so the correct statement is that 
g,h 6 SU (N c )/G c . The trick is now to find expressions for J± in terms of g, h which solve 
Equation (4.5). Indeed, it turns out that the desired expressions are [57] 

J+ = ^{g- 1 d + g-h- l d-h) (4.9) 
iN 

J- = -^(h- 1 d + h-g~ 1 d^g) (4.10) 

It is easiest to see these expressions satisfy Equation (4.5) in the axial gauge A- = 0, which 
can be accomplished by a gauge transformation that maps g —> gh^ = g' . Dropping the 
prime, one can easily check that J± do the job by using the identity 

D+ig^d+g) - d-ig^d+g) = (4.11) 

where D + is the adjoint covariant derivative. 

With these expressions in hand, we can rewrite Equation (4.4) as 

SW(g) = J d 2 x [J-6A+] = ^J d 2 xd-(g- l d + g)5gg- 1 . (4.12) 
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The remaining step is to find the action W(g) whose variation is Equation (4.12). As it 
happens, there is no manifestly SL^A^y-invariant 2D action for g which reproduces the 
expected SU(N c )a variation (i.e. the anomaly) which obeys Equation (4.12). However, 
it is possible to write an action with manifest SU(N c )y invariance and desired variation 
properties using the trick of allowing terms in the action to be written as integrals over a 
three-dimensional manifold which has the 2D spacetime as a boundary. The desired action 
is the Wess-Zumino-Witten action [3, 57] 



J d'xTrid^d^g) + ^ J^xdCe^Trig^dAgg^dBgg^dcg) 



(4.13) 



where the integral in the third term runs over D = IR 1 ' 1 x [0, 1], which in Euclidean space 
with ]R 2 S 2 would simply be a solid ball S 2 x [0,1]. The extension of the map g : 
S 2 — > SU(N C )/G C to S 2 x [0, 1] ->■ SU(N C )/G C does not have any topological obstructions 
since tt2(SU(N c )/G c ) = 0, and the topologically distinct extensions to D, classified by 
tts(SU (N c ) /G c ) = Z, all make the same contribution to the path integral provided N c £ Z. 
It is also straightforward to switch back to a general gauge as explained in [70], getting a 
result that depends on g and h and turns out to be writable as W(glv). For what follows 
we assume the switch to the general gauge has been made, rename gfv = g' and drop the 
prime. 

Finally, having obtained an expression for W(g), we can insert it into the path-integral 
bosonization/dualization machinery of [4-6, 69], which re-expresses Z{A) as a functional 
integral over g, so that 

Z(A)= J dgexpf^ J d 2 xTr{D ^g' 1 D» g) + i N c f(g, A)j , 

f (g, A) =-±- J d 2 x d£ e ABC Tr(g~ l d A gg- l d B gg- l d c g)+ (4.14) 

- i- J d 2 xe^Tr[iA^g~ 1 d u g -gd u g~ l + ig~ l A v g)\ 

where = + i[A^, •], and g £ SU(N c )/Zn c . This is the desired bosonized expression 
for adjoint fermions coupled to SU(N C ) gauge fields [72]. 

By construction, the bosonized theory transforms in the same way as the fermionic 
theory under local SU(N c )l x SU(N c )ji transformations, which act on g as g — > LgR -1 . 
Just like in the fermionic description, only the vectorial transformations (i.e. L = R) can 
be gauged. Of course, to see the non-invariance of W under local axial transformations 
requires a one-loop calculation in the fermionic language, while in the bosonic description 
this can be seen at the classical level. This is one of the key features of bosonization: 
one-loop physics in the fermionic language turn into classical physics in the bosonized 
description. 

Finally, we note that the fact that g lives in SU(N C )/G C rather than SU(N C ) is im- 
portant in ensuring that the bosonized action does not have a non-trivial symmetry under 
multiplication of g by elements of G c , which is fortunate since there is no such symmetry 
in the fermionic theory. 
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Transformation Properties 



Symmetry 




^R <P 


9 


SU(N C ) V 




Vip R V* V4>V* 


VgV^ 


SU(N C ) A 


Aij; L A^ 


Aty R A (j) 


AgA 


SU{N C ) L 




^R (f> 


Lg 


SU(N C ) R 




RiPrB) <p 


gR ] 








9 



Table 1. Transformation properties of matter fields in the bosonic and fermionic 2D theories. 



4.2 Mass and Yukawa terms 

To bosonize the theory described by Equation (2.5), we need to find a bosonization 'dic- 
tionary' for massive adjoint fermions coupled SU(N C ) gauge fields and scalar fields. Hence 
we now discuss the bosonization of mass and Yukawa-like terms in fermionic theories. 

The desired entries in the bosonization dictionary are most easily determined using 
spurion analysis 12 . The most general Yukawa-like term which can be made consistent with 

tp L -> Lip L L ] , 
ip R -> R^ R R ] , 

with L,R £ SU(N c )l x SU(N c ) r [the full symmetry group of Equation (4.7)] can be 
written as 

C F , y = Tr^A^B + h.c. = Tri> L A^ R B + h.c. (4.15) 

where A, B are spurion fields transforming as A — > LAR},B — > RBL^. There is also a 
symmetry A — > zA, B — > z _1 B, z G C\{0} which acts only on the spurion fields. Note that 
such a term is a relevant deformation of the fermionic action, since the mass dimension of 
A, B must be such that [A] + [B] = 1, and for particular choices of A, B one recovers the 
mass terms and (^-coupling terms in Equation (2.5). 

In the bosonized description of the theory the dynamical field is g, which transforms 
under SU(N c )l x SU(N c ) r as g — > LgR) . Since Cf, v is a relevant deformation of the 
fermionic theory, the same must be true for the term(s) induced by turning on the spurion 
fields in the bosonized theory. As it happens, the only relevant term involving g, A, and B 
which is allowed by the symmetries of g and the spurion fields is 



£-B,y = fJ- 



Ti(gB)TY{g ] A) + h.c] (4.16) 



where // is an undetermined positive parameter with mass dimension [p] = 1. The sign of 
fi is determined by the demanding that if A, B are set to constant values, as would be the 
case if one was bosonizing a mass term, so that A = B = y/m, m > 0, the bosonized theory 
must be stable, with a non-tachyonic mass term for fluctuations of g around the identity. 



12 For an interesting alternative prescription for the bosonization of mass terms of adjoint fermions see 
[73]. The relation between our results and those of [73] is unclear to us. 
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Bosonization Dictionary 



ja 


Tr if) , y + t a 'ip 


^Tr gid+gt" 


Jt 


Tr ip^f-t a il> 


^Tvd-g^gt a 
H [Tr(gB)Tr(g^ A) + h.c] 


^Yukawa 


Tr tpAipB + h.c. 



Table 2. Bosonization dictionary for SU(N C ) adjoint fermions with Yukawa and mass terms. 



The other relevant term consistent with the SU(N c )l x SU(N c )ji symmetry, which is 

(4.17) 

is ruled out by the U(l) symmetry acting on A, B. Finally, we note that the fact that the 
bosonized action is determined only up to a dimensionful parameter reflects the fact that 
the fermionic path integral with mass or Yukawa terms turned on is only really well-defined 
once a particular regularization scheme is chosen, leading to the fact e.g. {ipip) is scheme- 
dependent, and is a standard feature of bosonization both in Abelian and non-Abelian 
cases. The fact that this issue shows up already at the level of the classical action in the 
bosonized theory is another reflection of the fact that one-loop physics in the fermionic 
language becomes classical physics in the bosonic theory. 



B,y 



Tr{gAgB^) + h.c. 



5 Prom 3D gauge theory with fermions to a bosonized 2D theory 

With the results assembled in Section 2 and Section 4, we finally have all of the ingredients 
we need to bosonize 3D SU(N C ) gauge theory with adjoint fermions at large N c using the 
dictionary summarized in Table 2. Aside from the mass and Yukawa terms coupling ip with 
(p, the action of the bosonized description of Equation 2.5 will be given by Equation (4.14). 
Our task in Section 5.1 is to derive the terms in the bosonized action coupling g with (ft 
which aredue to the couplings between ip and <p in the 2D theory. Finally, in Section 5.2 
we put everything together to exhibit the full bosonized action. 

5.1 Couplings of <p and g 

To make use of the analysis carried out in the previous section we need to rewrite the 
fermionic interactions in the spurion formalism. For the once-and-future kinetic term, we 
have 

^Tr ([V^hsiV^}) 

Anticommuting the fermions, using the unitarity condition: (fxj)^ = 1 and the cyclic prop- 
erties of traces, we obtain the form: 

^Tr (^j L (p^ R (p - i> L cj)'ip R (f>^ , (5.1) 

and these two terms are separately Hermitian (recall that we are dealing with Majorana 
fermions). We can repeat the same procedure for the Wilson and mass term: 
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and using the same properties as above we get 

Tr (—^l^r ((> + —$ L (/> ipRtf ~ 2 ( — + m ) ) . (5.2) 



At this point we have to identify the spurion fields allowed by the transformation 
properties given in Table 1. In the generic case with both m,r/0we can introduce two 
spurions: one for the combined mass term containing the bare mass and Wilson mass and 
a second for the once and future kinetic term together with the Wilson term 13 

Tr (^ L A 1 ^ R B 1 + tp L A 2 ip R B 2 + h.c.) (5.3) 

with 

A x = ^am + rc B 1 = ^y/am + rc 

The fact that we need only two spurion fields, and not e.g. three, can also be seen by using 
the symmetries of the original fermionic theory such as ipi, — > —ipL, <\> ~ t 4> > m ~ * — m , r 
—r. 

Using the machinery built up in the previous section we arrive at the bosonized coun- 
terpart of the desired terms: 



8C = (m + -) W Tr 5 Tr 5 t _ ^ Tr(gt0)Tr(#t) _ M2c(r + 1) f ^ 

a la la 

(5.5) 



Had we started with the standard value of r in lattice calculations, r = 1, then we could 
have used just one spurion consistently with all the transformation of Table 1: 

Tr {^ L A^ R B + h.c.) (5.6) 

with 



yielding a particularly simple contribution to the bosonized Lagrangian: 

SC. = (?2±±) ^Tr 5 Tr fl t _ ^Tr^Tr^t) . 
\ a J a 



(5.8) 



Finally, we note that fj,,fii,fi 2 must scale as N® in the large iV c limit, due to the presence 
of two color traces in Equations (5.5) and (5.8), so that the contributions of these terms 
to the bosonized action have the same large N c scaling as the corresponding terms in the 
fermionic action. 



13 Note that as is explained in e.g. [39], the Wilson term coefficient should be chosen to obey |r| < 1 to 
ensure reflection positivity. 
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5.2 Bosonized form of 3D YM with an adjoint fermion 

We now have all of the ingredients we need to bosonize the theory described by Eq. (2.5). 
Making the simplifying choice r = 1, and combining the results above with those of Sec- 
tion 4.1, the bosonized version of S2D- 

s -=/ A { Tr (if F2 + ^f |0 ^ |2 + s |D " 912 ) 

+fh 2 TV g Tr g ] - -^Tr(#0)Tr (c/V) \ + NT(g, A) (5.9) 
a 1 J 

where we have defined rh 2 = Li(ma + c)a~ l and c = cfia. Together with the dimensional 
reduction formula in Eq. (2.5), this is our main result. The fact that bosonization is an 
exact procedure in two dimensions together with the large N c equivalence of Eq. (2.5) to 
the original 3D gauge theory implies that large N c 3D gauge theories coupled to adjoint 
fermions can be bosonized, with the resulting gauge theories written in terms of bosonic 
fields living in two dimensions. 

As the notation we chose above is meant to suggest, rh is a mass parameter which is 
related to the original bare fermion mass M, while c plays the same role as the 'speed of 
light' tuning parameter c in the dimensionally reduced theory. Indeed, as we emphasized 
in Section 2.1, the values of m,c (and hence rh,c) which correspond to a Lorentz-invariant 
continuum limit with a given physical quark mass must in general be determined non- 
perturbartively when the theory is strongly coupled, and get corrections order by order 
in perturbation theory when the theory is weakly coupled. This perspective makes the 
appearance of the parameter ix in the relations between rh, c and m, c especially natural, 
since the bosonized action encodes physics which is arises at the one-loop level in the 
original fermionic action. 

For the sake of completeness we also present the generic case with m/0,r/0: 

s = I Sx { n {wf + ^ D "* 2 + 5F M ) (5 ' 10) 

+m' 2 Tr 9 Tr 9 t - g( ^~ 1] T,(gU)1t(g^) - ^ TVfa^TrfeV )} + NT(g,A) 

where rh' 2 = Li\(ma + cr)a~ l and c' = CLt2<i, with the same comments concerning the 
connection of rh' , c' to m, c as above applying. 

6 Conclusions 

The recent realization that D-dimensional SU(N C ) gauge theories compactified on spatial 
tori e.g. R 1 ' 1 x S D ~ 2 with adjoint fermions stay in a center-symmetric phase for arbi- 
trary torus size L has important consequences for both analytic and numerical studies of 
gauge theories. Analytic developments have so far mainly focused on the LA^A st rong <C 1 
regime, which is center-symmetric but not volume-independent, and have very recently 
lead to dramatic progress toward a non-perturbative continuum definition of a class of 
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gauge theories [74, 75]. The main applications of the notion of large N volume indepen- 
dence, which holds when LNA strong ^ 1> has so far been for Monte Carlo simulations, 
since it means Eguchi-Kawai reduction to one plaquette is valid for such theories [76-84]. 
Our results here suggest that the notions of large N c volume independence and large iV c 
dimensional reduction may also be quite useful for obtaining novel analytical results. 

Our main results are Equation (2.5) and Equation (5.9). The notion of large N c 
dimensional reduction has appeared before in [37, 56] in the context of reduction from a 
4D theory to a 3D one, but we believe the presentation here clarifies a number of aspects 
of large N c dimensional reduction. As an application of this technology, we showed that 
(perhaps unexpectedly) 3D SU(N C ) gauge theory with one flavor of adjoint fermions can 
be bosonized, with the bosonized action given in Equation (5.9). A remarkable feature 
of this result is that the bosonized action corresponding to the large N three-dimensional 
theory is local, but lives in two dimensions. 

There are a large number of possible applications and extensions of our results. For 
instance, while we focused our discussion on 3D SU(N C ) gauge theory with Nf = 1 adjoint 
fermions, the techniques we discussed are clearly more general. For instance, one could 
bosonize 4D gauge theories by working on R 1 - 1 x S 1 x S 1 , discretizing on a lattice with V 2 
sites, performing the large N c dimensional reduction, and then applying a generalization 
of the adjoint fermion bosonization dictionary to multiple flavors to find the 2D bosonic 
actions corresponding to 4D large N c gauge theories. It would be very interesting to work 
out the explicit form of such bosonized theories. 

A crucial issue which deserves further studies is the analysis of the spectrum beyond 
the quadratic level. Clearly this is a difficult task but in principle we could combine 
the expansion around the center symmetric vacuum with the powerful tool of discretized 
light-cone quantization (DLCQ) to analyse the full interacting spectrum. The idea is to 
compactify also the light cone space x~ so that all the momenta will be quantized and 
the 2D Hamiltonian could be diagonalized numerically at large- N using only single trace 
states. Previous studies (see e.g. [85-89]) were always carried out in the trivial vacuum, 
while in our analysis of the dispersion relations in the free case, we saw that the presence of 
a non-trivial Polyakov loop turns out to be crucial to encode the momenta in the reduced 
direction. Hence it would be especially interesting to repeat the DLCQ calculations using 
the center-symmetric vacuum. 

Another interesting direction is to start with D-dimensional gauge theories with ex- 
tended supersymmetry, and attempt to define lattice regularizations on T D ~ 2 which pre- 
serve some supersymmetry at finite lattice spacing (see [90] for a review of lattice super- 
symmetry) . This should result in a version of large N c dimensional reduction with manifest 
supersymmetry and might enable some interesting calculations which would be analytically 
intractable using the simple-minded regularization we adopted here. It should also be very 
interesting take advantage of the super-renormalizability of 3D gauge theories to do the 
perturbation theory calculations necessary to find the expressions for m, c which lead to 
the desired continuum limit [47-53]. 

As a final example of an application of our results, one could try to extract some 
phenomenology from the bosonized form of the 3D theory, or the appropriate 4D general- 
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izations. For instance, it is possible to use the methods of [7, 91-93] to compute the fc-string 
tension of the 3D theory in the L ~ iV" 1 large N c limit, as we will discuss elsewhere [94]. 
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